Band structure and gaps of triangular graphene superlattices 
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General properties of long wavelength triangular graphene superlattice are studied. It is shown 
that Dirac points with and without gaps can arise at a number of high symmetry points of the 
Brillouin Zone. The existence of gaps can lead to insulating behavior at commensurate fillings. 
Strain and magnetic superlattices are also discussed. 



I. INTRODUCTION 



Graphene is a two dimensional metal when carriers are 
induced by an electric fielctirii. A gap at the Fermi level 
has been observed by STM measurements^!^ (see also^). 
We analyze the gaps induced by a periodic structure, 
and the possibility that these gaps are generated sponta- 
neously. 

Graphene superlattices have been observed in 
graphene layers grown on transition metals^i^ (see 
alsoiS^J^). Superlattices are also found in graphene 
grown by the decomposition of SiGii. In general, 
graphene superlattices can have interesting properties, 
such as highly anisotropic transport properties'^, or 
Dirac points at finite energiesi^"— . In general, the study 
of the properties of graphene superlattices has attracted 
great interest, due to the many novel features their elec- 
tronic spectra can show^ir— . In the following, we an- 
alyze general properties of the spectra of graphene su- 
perlattices with a two dimensional triangular periodic- 
ity. These supperlattices share the symmetries of the 
graphene lattice, and are commonly found in graphene 
layers grown on metallic surfaces. As discussed below, 
these superlattices can show a gap at the Fermi energy 
for a number of commensurate fillings. It seems likely 
that they can be formed spontaneously on very uniform 
substrates, such as BNi, or due to intrinsic instabilities 
of graphene. 

In this paper, we study the general properties of tri- 
angular graphene superlattices created by a scalar po- 
tential, followed by a discussion of strain and magnetic 
superlattices. 



II. BRILLOUIN ZONE OF TRIANGULAR 
GRAPHENE SUPERLATTICES 



We define the lattice vectors of the graphene lattice as: 



ai 



0-2 = -j^x 




FIG. 1: (Color online). Examples of Brillouin Zones of super- 
lattices. Left: ni = l,n2 = 4. Right: ni =2,n2= 4. 



A triangular superlattice is described by the unit vectors: 



hi = mai -I- 1120,2 

1)2 = —112(11 + (ni + 7l2)a2 



(2) 



where ni and n2 arc arbitrary integers different from 
zero. 

There are three types of high symmetry points in the 
Brillouin Zone of a triangular lattice, F, M and K. There 
are two inequivalent K points, K and K' , at the corners 
of the hexagonal Brillouin Zone, and three inequivalent 
M points, at the centers of the edges. Time reversal 
exchanges K and K' , while leaving the F point and the 
three M points unchanged^. The vectors which define 
these points are such that: 
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The low energy states of graphene lie close to the K 
and K' of the original Brillouin Zone. The positions of 
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these points in the supcrlattice Brillouin Zone is deter- 
mined by: 



2tt -> ^ A-K 
—712 Ksh2 = ^"^"'^a 



Ksh = —ni 
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27r 
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Thus, when 27ii + 712 is a multiple of three the graphene 
K and K' points will be mapped onto the Ts point of 
the superlattice Brillouin Zone. Otherwise, they will be 
mapped onto the corners of the Brillouin zone, Ks and 
K'g. Examples of superlattice Brillouin Zones are given 
in Fig. m 



FIG. 2: (Color online). Points in the superlattice Brillouin 
Zone where degeneracies occur when Vc 7^ 0. 



distances from the points, we can write an effective 
hamiltonian: 



III. DISPERSION NEAR HIGH SYMMETRY 
POINTS 



/ vfG_ 



(7) 



The model 



which gives an anisotropic Dirac equation with a gap: 



We study superlattices induced by a modulation of the 
on site energy of the tt orbitals. We assume that it is a 
weak perturbation of the graphene Dirac equation, ex- 
cept in cases where degeneracies occur. We consider the 
Fourier components of the potential with lowest wavevec- 
tor, of modulus: 
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We write the potential as the sum of s symmetric part Vg, 
and an antisymmetric part, Aq, with respect to the inter- 
change of sublattices. We neglect for the moment short 
wavelength components which mix the two inequivalent 
Dirac points of the unperturbed graphene layer. 

We analyze the changes in the Fermi velocity near the 
Dirac energy induced by the superlattice potential, and 
the points in the lowest bands of the superlattice where 
degeneracies persist when Vg ^ and Aq = 0. As dis- 
cussed below, this situation gives rise to a new set of 
Dirac equations at finite energies. 
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At the Ks and K'g points, there are three degenerate 
levels for Vg = Aq = 0, with energy = vpG/\/3. 
When Vg 7^ these three levels are split into a doublet, 
with energy ej^^ ~ vpG/^/S + Vg/2, and a singlet, at 

e|fg = vpG/y/S — Vg- Expanding around the Ks point, 
the effective hamiltonian for the doublet is: 



( VfG , Vg _ vpk^ vphy _ 3iAc 

\/3 ^ 2 2 2 4 
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This is the two dimensional Dirac equation with a mass 
term. The dispersion relation is: 
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B. Dirac energy at the Ts point 

We consider the case when the K and K' points of the 
graphene Brillouin Zone are mapped onto the Ts point 
of the superlattice Brillouin Zone. Using lowest order 
perturbation theory, we find, near the Ts point a renor- 
malization of the Fermi velocit y ^^1-'^^ : 
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There is a twofold degeneracy at the three Ms points, if 
Ag = 0. The energy of these states is vpG/2. At finite 



There are two sets of degenerate bands, derived from the 
K and K' points of the Brillouin Zone of graphene. This 
degeneracy will be broken by short wavelength terms in 
the superlattice potential. 



C. Dirac energy at the Ks and K's points 

The renormalization of the Fermi velocity near the Ks 
and Kg points is the same as in eq. H) 

There are doubly degenerate states, even when Vg ^ 0, 
at the six inequivalent points at positions Lg = Ks/2, 
as shown in Fig. The energy of these points is 
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FIG. 3: (Color online). Top: Low energy bands for a 12 x 12 
superlattice, with Vg ~ 0.3eV and Ac = 0. Bottom: Left: 
Details near the Ms point. Right: Detail near the Ks point. 
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FIG. 4: (Color online). Top: Low energy bands for a 12 x 12 
superlattice, with Vg ~ 0.3eV and Ac = O.leV. Bottom: 
Left: Details near the Ms point. Right: Detail near the Ks 
point. 



= vfG/2. Expanding around these points, we find 
an effective anisotropic Dirac equation, given by eq.(I7]). 

There are another set of doubly degenerate states at 
the Ms points. The two states arise from the K and 
K' points of the original graphene Brillouin Zone. The 
degeneracy persists when Vq ^ and Aq ^ 0, and is 
only broken by short wavelength components of the su- 
perlattice potential. When these components are finite, 
an effective anisotropic Dirac equation will arise similar 
to that in eq.d?]). 

For Vg = = there are six degenerate states at 
the Ts point. The long range part of the superlattice 
potential will hybridize states which are derived from the 
K and K' points of the original graphene Brillouin Zone. 
We obtain two sets of isotropic Dirac equations, described 
by eq.®, and two degenerate states. The short range 
part of the Dirac equation will break these degeneracies. 
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FIG. 5: (Color online). Density of states for a 12 x 12 super- 
lattice. Red: Vg = 0.3eV and Ac = 0. Blue: Vg = 0.3eV 
and Ag = O.leV. 



D. Results 

We analyze the bands induced by a iV x superlattice. 
The hopping matrix between tt orbitals in neighboring 
carbon atoms is t — 3eV. The bands for Vg — 0.3cV and 
Ag = are shown in Fig.|3l The bands show Dirac points 
at the Ms and Ks points. When Ag is increased to 
Ag = O.leV a gap appears between successive bands, as 
shown in Fig.|31 The density of states for those two cases 
is shown in Fig. [5] Note that the potential breaks the 
electron-hole symmetry of clean graphene, and the gaps 
are not of the same magnitude for positive and negative 
energies. 



The results are in reasonable agreement with the ana- 
lytical description in the previous section. A gap of order 
2Ag is induced at the Ms point. In order for this gap to 
be possible, the following inequalities must be satisfied: 
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The scaling properties of the Dirac equation imply that, 
if the dimension of the superlattice is increased, G — >■ 
AG, with A < 1, a rescaling of the superlattice potential, 
Vg — AVg,Ag —J- AAg, will lead to the same band 
structure, with energies scaled as, E ^ \E. 
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IV. STRAIN SUPERLATTICES 

A superlattice can also be produced by inducing 
strains, which modulate the interatomic hoppings. The 
corresponding perturbation can be seen as a gauge field, 
A, which shifts locally the momentuni^. A simple case 



is when the strains are due to height modulations, h{r), 
which can be induced by a substrate. In terms of the 
Fourier components of the modulation, h^, and allow- 
ing for the relaxation of the in plane displacements, the 
effective gauge field can be written as^: 
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(12) 



where the tensor h^i is the Fourier transform of the func- 
tions: 



iG, 



dh dh 

dxi dxi 



(13) 



and A and /i are the elastic Lame coefficients of graphenc. 
The field in eq. [T^] has opposite signs in the two valleys 
in the Brillouin Zone. 

The calculation of the effective magnetic field induced 
by the gauge field in eq. [T^ is simplified when, as 
in the previous sections, only one component, hir) ~ 



e'"'"", in a superlattice is considered. The 



^GEi=l....,6 

tensor in eq.[l3]has non zero components for all combina- 
tions of the type Gfe + G;. When the gauge field, A(G) is 
parallel to G), the vector potential can be gauged away, 
and does not induce an effective magnetic field. This 
implies that out of the 18 possible values of the vector 
Gfc + Gj, only six contribute to the effective magnetic 
field. These vectors are given by G = G(3/2, 73/2) and 
the vectors equivalent to it by a symmetry transforma- 
tion. After some algebra, we obtain for the effective mag- 
netic field: 
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.(f) 



^^^27a/3^3^2 
a A + 2^ 8 



G 



2 cos( VSGy) -I- 4 cos 




(14) 



where /3 = 91og(i)/91og(a) « 2 — 3, t « 3eV is the near- 
est neighbor hopping, and a ~ 1.4A is the distance be- 
tween nearest neighbor carbon atoms. The superlattice 
defined by the effective magnetic field has a unit vector 
of length -\/3G, so that the area of its unit cell is a smaller 
than the area of the unit cell of the original superlattice 
by a factor 1/3. 

Using eq. HH and G = 27r/L, where L is the length of 
the unit vector of the superlattice, we can write the mag- 
netic length associated to the maximum effective field in 
the system as: 



1 A + /i /3 9V3tt^ hi 



A + 2/i a 



L3 



(15) 



where h^ax — 6/1^ is the maximum value of h{G), as- 
suming Hq > 0. For values h^ax ~ Inm, and L « 40nm, 
we find Ib ~ 14nm, so that the effective field is such that 



Tjmax ^ 1 
strain ^ 



2T. 



MAGNETIC SUPERLATTICES 



A superlattice can also be induced by a spatially mod- 
ulated magnetic field. A combination of a modulated 
magnetic field and a scalar potential opens a gap at the 
Dirac energy, and the resulting insulator is a Quantum 
Hall systen*2i^ with chiral currents at the boundaries^. 
Here we obtain this effect using second order perturba- 
tion theory, instead of the arguments used iii^l. As in 
the previous sections, we assume the simplest periodicity 
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compatible with the superlattice hexagonal symmetry 



V{v)=Vg J2 

1=1, ■■■ ,6 

B{v) = Bg 

1=1, ■■■ ,6 



(16) 



The cigcnstates of the unperturbed hamiltonian at the 
K and K' points of the Brillouin Zone can be written 
as |0)a = (1,0) and |0)b = (0,1), which each compo- 
nent of the spinor corresponds to one sublattice. These 
states arc hybridized with states \G±)k = (l,±e"^e) 
and \G±)k' = (1, Te^*'^^), with energies eg = ±uf|G|, 
and e^-^G = {Ga^ + iGy)/\G\. 

The energy of states \0)a and |0) b are modified in dif- 
ferent ways by virtual hoppings into states \G±)k and 
\G±)k', leading to gaps in both valleys. Moreover, the 
gaps have different signs, showing that time reversal sym- 
metry in the system is broken, and that a Quantum Hall 
phase has been induced. The gap can be written as 



A = ± E 



1=1, ■■■ ,6 





A^{Gi) + iAy{Gi) 











(17) 



orbitals. An electron in a given carbon atom experiences 
the image potential due to the screening. A change of 
position of Az leads to a change of the image potential 
of order ee* Az / {Ad'^) , where e* = e(eo - l)/(eo + 1) is 
the image charge, eo is the dielectric constant of the sub- 
strate, and d is the distance to the substrate. A vertical 
displacement of Az ~ lA when the graphene layer is at 
a distance d ~ SA of the substrate can lead to shifts of 
the onsite energies of order O.leV. The electronic gain of 
energy due to the existence of a gap, per unit cell, is then 
Eeiec « Ag « ee*Azd-2. 

The elastic energy per unit cell required to create a 
staggered distortion of amplitude Az is of order Eeias ~ 
kAz^o^^, where k sa leV is the bending rigidity of 
graphene. 

A gap will exist above a threshold for the superlattice 
potential, Vg ^ Aq > vfG{1/V3 - 1/2). The area of 
the Brillouin Zone of the supcrcell, \/3G'^/2 should be 
close to the area within the Fermi surface of the unper- 
turbed graphene, irkp. Hence, a spontaneous staggered 
distortion is favored if: 



E, 



elec 



where A^^ is the vector potential, which we define as: 



f Eelas 
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ee 
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d^ 
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(20) 
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(18) 



where $o — dnj c is the quantum unit of flux. Using this 
expression, we finally obtain: 



A = 12 



BgVg 



(19) 



in agreement with21. 



VI. SELF CONSISTENT OPENING OF A GAP 

The previous analysis shows that gap can open at finite 
energies in graphene in the presence of a superlattice po- 
tential with a staggered component. When the number 
of carriers is such that only a small number of subbands 
are completely filled and the rest are completely empty 
the electronic energy will be lowered in the presence of 
the gap. A lattice distortion which leads to the appro- 
priate potential will be energetically favorable if the gain 
in electronic energy exceeds the formation energy of the 
distortion, as in the Peierls instability in one dimension. 

In graphene on top of a metal or other substrate with 
a large dielectric constant, as in^, out of plane displace- 
ments lead to changes in the on site energies of the tt 



The last equation in (j20|) implies that G oc kp, and from 
the first two we obtain that vpG < {ee*)'^a'^d^'^K^^ . 
Hence, a spontaneous distortion is energetically favored 
for carrier densities such that: 



^'^F ^ „2 ,.2 J8 



(21) 



and the staggered distortion, with inverse wavelength 
G ^ kp oi the order: 



Az 



vpkpd'^ 

ee* 



(22) 



For a metal, we have eo = oo and e* — e. For 
a sa 2A, d w sA and n ~ leV, we find that a staggered 
corrugation is energetically favored for carrier densities 
P ^ lO^'^cm^^^ leading to maximum deformations of 0.2 
A. For Si02, where eo = 3.9 and e* « 0.6e, the corruga- 
tions takes place for carrier densities /o < 3 x lO^^cm"^^ 
and maximum deformations of 0.15 A. 

The formation of these long wavelength modulations 
is only possible in systems with a high degree of order, 
as the superlattice features will be reduced by disorder. 
For instance, weak scatterers with concentration iiimp, 
which change the onsite potential by an amount of order 
Ag, lead to a mean free path /gj ~ Vp / (A^kpa'^nimp) ■ 
A large concentration of these defects, n 



imp 



will 

suppress the formation of a superlattice, while affecting 
little the conductivity, a ^ (e^/h) x (kplei) ^ (e^/h) x 
iv'pa-'A^'). 
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VII. CONCLUSIONS 
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